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Abstract. We consider the running of the spectral index as a probe of both inflation 

itself, and of the overall evolution of the very early universe. Surveying a collection 

w^ ■ of simple single field inflationary models, we confirm that the magnitude of the 

running is relatively consistent, unlike the tensor amplitude, which varies by orders 



O ' of magnitude. Given this target, we confirm that the running is potentially detectable 

■j-j , by future large scale structure or 21 cm observations, but that only the most futuristic 

jrt I measurements can distinguish between these models on the basis of their running. For 

any specified inflationary scenario, the combination of the running index and unknown 
/»f-N , post-inflationary expansion history induces a theoretical uncertainty in the predicted 

^ ' value of the spectral index. This effect can easily dominate the statistical uncertainty 

Qy . with which Planck and its successors are expected to measure the spectral index. More 

positively, upcoming cosmological experiments thus provide an intriguing probe of 
^f) • physics between TeV and GUT scales by constraining the reheating history associated 

I — , with any specified inflationary model, opening a window into the "primordial dark 

f^ ' age" that follows the end of inflation. 
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1. Introduction 

Concordance cosmology fits the overall properties of the universe in terms of a handful of 
physical parameters. This set is not fixed, but is the minimal number required to match 
the data [UEIE]. The concordance parameter set will expand as more subtle physical 
effects are resolved in the data, while current variables may someday be determined 
independently of cosmological observations. Within the concordance model we may 
delineate sectors, or subsets of related parameters. For example, Qb, ^c and ^a (the 
relative contributions of baryons, cold dark matter and dark energy to the present- 
day density of the universe) define the composition of the universe. These quantities 
stand in for our ignorance of aspects of fundamental physics: fib, ^c and fl\ are free 
parameters because we do not understand baryogenesis, the dark matter abundance, 
and the amplitude of the vacuum energy. 

Two concordance parameters specify the primordial perturbation spectrum - the 
amplitude Ag and spectral index Ug. If perturbations are generated during inflation, Ag 
and Us form the nucleus of an inflationary sector. Many further observables are related 
to the primordial perturbations, including tensor modes (their amplitude relative to the 
scalar perturbations r, and their spectral dependence, nt), non-Gaussianity, running in 
the scalar spectral index (a^ = dris/dlogk, where k is comoving wavenumber), features 
in the power spectrum, and relics generated at the end of the inflationary epoch. These 
"optional" observables are fingerprints of specific inflationary scenarios, in that most 
models predict that most of these parameters are very small: finding a non-zero value 
for any one of these quantities would slash the number of viable inflationary models. 

Predictions for ra^, a^, r, et al. define a mapping from the fundamental physical 
description of inflation (namely, the gravitational action and relevant field content, 
which is the inflationary sector of very high energy particle physics) into the parameter 
space of observables. The canonical example of this mapping is the "zoo-plot" which 
locates single field, slow roll inflationary models in the {ns,r) plane [1]. Unfortunately, 
the clarity of this plot is partially due to its mapping a subset of inflationary models 
into a subset of inflationary observables. Furthermore, as we explore below, our lack 
of knowledge of the post-inflationary expansion history renders this mapping inherently 
imprecise. 

The first WMAP data release marked a turning point in constraints on inflationary 
scenarios [5]. A host of currently ongoing and planned observations, looking at both the 
CMB and complementary datasets, ensures that the de facto observational campaign 
to constrain the inflationary era will last for at least the next several decades. The 
purpose of this paper is to highlight the importance of the running, or scale dependence 
of the spectral index, as a key inflationary parameter. Our first task is to explore the 
correlation between ag and the measured value of Ug induced by the (unknown) post- 
inflationary expansion history: the equation of state for the post-inflationary universe 
determines how rapidly modes reenter the horizon, and thus the moment during inflation 
when a given mode leaves the horizon. By definition, a^ fixes the scale dependence of 
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Ug, so the observed value of n^ depends on a^ via the post-inflationary expansion history. 
Consequently, assumptions about the post-inflationary expansion history of the universe 
add a (sometimes implicit) theoretical prior to constraints on inflationary parameters. 

We consider several explicit inflationary models, namely natural inflation (a 
sinusoidal potential) [6], 0" inflation, where n is a continuous parameter (motivated by 
models for which n takes fractional values [3 [8]), and "inverted" potentials of the general 
form A^ — XM^^^(j)'^/n [9[ [101 El]- This is not an exhaustive survey, but it captures 
many well-known models. We consider the theoretical uncertainty in predictions for 
inflationary observables via the matching equation, which connects comoving scales 
in the present universe to the moment during inflation at which the corresponding 
perturbations were produced. The matching depends on the rate at which modes 
reenter the horizon, which is in turn a function of the post-inflationary expansion history. 
However, there are few constraints on the detailed expansion history of the universe prior 
to nucleosynthesis, and the details of this epoch substantially modify the predictions of 
explicit inflationary models. 

We consider three specific scenarios for the post-inflationary expansion. The first 
is instant reheating, for which the inflationary observables are well defined. Secondly, 
we allow a matter dominated phase of unknown duration, followed by thermalization. 
Finally, we consider a post-inflationary equation of state fixed only by general physical 
arguments and direct observational constraints. The unknown expansion history leads 
to an uncertainty, AA^, in the number of e-folds, N. For simple inflationary scenarios, 
10~^ ~ \'^s\ ^ 10"'^, providing a benchmark value for this parameter. For a given AA^, 
the measured value of Ug changes by Aug ~ a^AA^. Assuming a long period of matter 
domination before thermalization yields AA^ ~ 10 and Aug ~ 0.005 (similar to the 
statistical errors expected from Planck), and AA^ and Aug can be much larger if we 
drop all theoretical priors on the post-inflationary dynamics. Note that this uncertainty 
is directly induced by the running: by definition, if a^ = 0, Ug is not a function of k, 
and its value at the pivot scale is unaffected by the details of reheating. 

Despite the large uncertainty in Ug induced by the interaction between the running 
and post-inflationary expansion, Og is undetectable by a cosmic variance limited CMB 
mission for simple inflationary models. However, the next decade will see dramatic 
improvement in the spatial volume and redshift range probed by large scale structure 
observations, and the beginning of the exploration of the dark ages via high redshift 
21 cm emission. These measurements are directly sensitive to the primordial spectrum 
over a broad range of scales, and can improve limits on Og. Unlike the CMB, probes 
of large scale structure sample the universe volumetrically, and are not subject to 
Silk damping which washes out the CMB at small scales. We explore the sensitivity 
of galaxy surveys and 21 cm experiments to a^, as a function of their basic design 
parameters, confirming that ambitious surveys can detect a^ at the level predicted 
by simple inflationary models. However, even highly ambitious proposals would 
have trouble measuring ag with enough precision to distinguish otherwise degenerate 
inflationary scenarios. On the other hand, the correlation between Ug and the post- 
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inflationary expansion history induced by the running means that measurements of 
Ug yield information about the post-inflationary expansion history following a specified 
inflationary scenario. This effectively uses the theoretical uncertainty in Ug to gain 
information about the post-inflationary equation of state for a given inflationary model. 
While this is a single number that depends on the integrated expansion history of the 
universe, it probes physical processes that occur between the TeV an GUT scales, which 
are otherwise almost unconstrained. Specifically, this number yields information about 
the mechanism of post-inflationary reheating, and whether a specific inflationary model 
is compatible with - or requires - a transient matter dominated phase in the early 
universe, as well as more exotic scenarios such as kination [12] (or deflation [13]) or a 
phase of cosmic string dominated evolution [T3] . 

Given that we are considering simple, single field models, our analysis touches on 
a great deal of earlier work. Large, negative values of Og signal either the imminent 
breakdown of slowroU, or that higher order terms in the potential dominate the 
inflationary evolution [151 HBl [E]- Conversely, a positive running can be associated 
with primordial black hole production [181 [H]) leading to further observational 
constraints. The correlation between the post-inflationary thermal history and 
inflationary observables is usually expressed via the matching equation, as surveyed 
in [20l [m [22]. Our analysis overlaps with discussions of the "likely" values of 
inflationary parameters: in particular, the natural value of r in inflation is hotly debated 
[231 [211 [251 [23 while Kinney and Riotto have discussed the theoretical uncertainty in Ug 
arising from higher order terms in the inflationary potential [27] . Our treatment of the 
joint constraints on ag from CMB, 21cm and surveys is complementary to treatments 
of combinations of CMB and gravitational wave experiments [28]. Further the ability 
of future 21 cm experiments to detect a running index was analyzed by [29], while [30] 
discusses the detectability of a running index through its effects on structure formation. 

The paper is organized as follows. In section [2] we review canonical single field 
inflation. We describe the matching of physical scales to field values during inflation 
and introduce the effective equation of state parameter w [311|32]. In section[3]we collect 
predictions for the perturbation spectra of three general classes of single field models: 
single field, natural, and hilltop or inflection point inflation. We explore the parameter 
uncertainty and model degeneracy induced by the unknown post inflationary thermal 
history. In section [H we explore the constraints that may be placed on the tilt and the 
running of the scalar spectrum by combinations of future astrophysical observations. We 
focus on observations of the cosmic microwave background, galaxy surveys and cosmic 
tomography using the 21 cm line of neutral hydrogen. In section [5] we employ the Fisher 
matrix formalism to forecast constraints on the post inflationary parameters for a given 
inflationary model. Finally, we conclude in section [61 
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2. Inflation, The Matching Equation, and the Primordial Power Spectrum 

2.1. Equations of Motion and Slow Roll Expansion 

We begin by gathering results for the dynamics and observable properties of single field 
inflation. In a flat Friedmann- Robert son- Walker universe dominated by a scalar inflaton 
with potential V{(f)), the Einstein and Klein-Gordon equations are 
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Variables have their conventional deflnitions and we use the reduced Planck mass, 
suppressing factors of vr throughout our analysis. 

We work in the potential slow roll (PSR) expansion - this is less elegant than the 
Hubble slow roll (HSR) expansion [33], but expresses physical observables directly in 
terms of the potential. The flrst three slow roll parameters are 
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2.2. The Perturbations 

Quantum fluctuations of the metric and matter flelds during inflation give rise to both 
curvature and density perturbations. These fluctuations are characterized by the gauge 
invariant variable (, which measures the density fluctuations on hyper-surfaces with 
uniform curvature. The power spectrum of ( can be solved in the adiabatic regime 
{e,r] ~ const.) as an expansion about the exact power law solution [3l] 
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where C = —2 + log 2 + 7 ^ —0.73 and 7 is the Euler-Mascheroni constant. Similarly, 
the primordial spectrum of gravitational waves is 

V,{k)= [l-{C + l)eHrT7T{i^] ' (7) 
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The spectra are conventionally written as power laws, 



k 



ns{k) — l 



Vdk) = A,{h) ( ^ ) , (9) 
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nt{k) 



Vu{k) = Mk,) [jj , (10) 

where Ug and Ut are the spectral indices. With a specified pivot A;^, we measure ^^(A;^), 

as = dns{k) /d log k\k=k^, and r(A;^) = Vh{k-^)/'P(^{k^,). For single field models, we have 
the consistency relation Ut = — r/8 (to first order in slow roll), and 

^(n,-l)= -3e + r]-(^ + Uc\e' + {8C-l)er] + ^r]' 

C^-0e + ---, (11) 

'^'"^ -16er]-24e^-2^ + ---, (12) 



dlog/c 



r = 16e 



1 - ^e + ir/ + 2C(2e - r]) 



+ ■••, (13) 



where '■ ■ ■' denotes higher-order terms in the slow-roll expansion. 

2.3. Matching Equation, Thermalization, and Expansion History 

The spectral indices and running are measured at the pivot, k^,, a specific, physical 
scale in the present-day universe. To compare inflationary predictions to observed 
inhomogeneities, we match this length to the field value at which the expressions of 
the previous subsection are evaluated. During infiation, a mode leaves the horizon 
when its wavelength matches the Hubble radius: it reenters the horizon when the post- 
infiationary Hubble radius expands to this (comoving) scale. Comparing a comoving 
scale, k, with the current Hubble scale we solve for N{k), the number of e-folds before 
the end of infiation at which the mode left the horizon: 

.,.) ^10. (£-).. o.(^)-.o,(-A.). ,,, 

Assuming slow roU, we have the usual result 

N{k) ^ log (^) = T" H{t)dt ^ -^ /"' :^d0 . (15) 

For a specific infiationary potential, the ratio Hk/Hend can be accurately evaluated but 
N{k) also depends on the growth of the horizon, {aH)~^, from the end of infiation to 
the present day. We assume that at any given instant the post-infiationary universe has 
a well-defined equation of state, w, the ratio between pressure and density|§| 

§ This can be an effective parameter. During coherent oscillations, w rapidly oscillates between ±1, 
but the time-averaged value is w = in a quadratic potential. The growth of perturbations depends 
on the (time-dependent) ratio of the their frequency to that of the field oscillations [35], but our focus 
here is the background dynamics. 
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Nucleosynthesis (e.g. |i36j) and evidence for a cosmological neutrino background 
[37] imply that the universe was in thermal equilibrium at MeV temperatures, and 
the evolution of the universe between this epoch and the present day is assumed to 
be well-understood. The current microwave background temperature is C(10~'^) eV, 
from which we infer that the universe has expanded by a factor ~ 10^° since neutrino 
decoupling. However there is no concrete evidence that the universe was thermalized at 
higher densities. Conversely, as inflation ends, the universe contains few particles and 
is far from thermal equilibrium. The mechanism of reheating is unknown, and usually 
unspecified by the inflationary potential on its own. We parameterize our ignorance by 
specifying that the universe thermalizes at a temperature Treh ~ Prih- ^^ noted above, 
the only firm constraint on this parameter is pj^f^ > (9(10) MeV. If inflation is a GUT 
scale phenomenon, the epoch between the end of inflation and the MeV scale spans a 
range of 10^* in energy. Consequently, the ratio of scales between the end of inflation and 
neutrino decoupling can far exceed the separation between neutrino decoupling and the 
present day. Moreover, over most of this epoch the basic properties of particle physics 
itself are poorly understood, and amount to a "primordial dark age" . 

The "default" assumption is that the post-inflationary universe undergoes an 
effectively matter dominated period of coherent oscillations, followed by thermalization, 
followed eventually by matter-radiation equality and the onset of structure formation. 
However, a huge range of mechanisms could interrupt this simple picture, and their 
impact on the expansion rate can typically be described by the equation of state w, 
the ratio of the pressure to the density. We are thus implicitly limiting ourselves to 
a barotropic fluid, but the discussion here could easily be generalized. These range 
from kination {w = 1) ^2\, to frustrated cosmic string networks [H] {w = —1/3), or 
even a short burst of thermal inflation [38] during which some modes would be pushed 
outside the horizon for a second time. Consequently, we break the post inflationary 
evolution into two parts, the unconstrained post-inflationary era, and the hot big bang 
era (which runs through to the present day), during which the universe is thermalized 
and populated with familiar standard model particles. Considering the former, we write 



faendHend\ A\og{aH) ^ 

log Tl— = xn A log a. 

V arehH-reh J A log O 



(16) 



Rather than adding terms to equation ( TT6l) to account for each possible phase in the 
post-inflationary universe, we define an effective equation of state 



1 r 

w = -— / u;(a)dloga. (17) 

Introducing w replicates the growth of the scale factor and horizon during these 
intermediate stages, as sketched in Figure [D Using standard results for the evolution of 
the horizon and scale factor with a fixed equation of state we derive 

(aendHend\ {I + ?>w) f 2 Preh\ ,,oN 

log ^r~ = ~7^7^~^ log Q T^ • ^^ 

\arehHrehJ 6{1 + W) \SVendJ 

The logarithmic term follows from noting that a = p + 3p = at the end of inflation, 
giving Pend = 3Knd/2. 
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Figure 1. The evolution of tiie Hubble horizon for the inflationary universe is 
shown in cartoon form. The x-axis describes the cosmological scale factor a{t) on 
an approximately logarithmic scale. The parameter w describes the growth during the 
epoch between the end of inflation with the standard hot big bang era. Note that for 
smaller w the value of N at which the pivot leaves the horizon is similarly decreased. 



Almost any conceivable evolution between the end of inflation and the present day 
is now encoded in w and preh- This argument relies on astrophysically relevant modes 
being outside of the horizon until after thermalization: modes inside the horizon respond 
directly to the equation of state, and need the full transfer function - this distinction 
will be important for direct detection gravitational wave experiments [391 SO] • Assuming 
reheating is followed by the radiation and matter dominated epochs, one arrives at the 
matching equation: 



N{k) =56.12 -log 




(19) 



where k^ = 0.05 Mpc~^ is chosen as the pivot scale in what foUowsl]]] Notice that 
all dependence on the present-day Hubble parameter, h, cancels. A fluid with p > p 
has a superluminal sound speed, so we expect w < 1. Conversely, inflation ends at 
Vend, SO w > —1/3, to cusurc that the comoving horizon does not shrink further during 



Note that the optimal choice of pivot depends on the dataset being considered [U] 25] . 
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Figure 2. We illustrate the impact of changing our the assumed value of Preh- 
The left hand panel demonstrates that setting a high value of preh reduces the total 
uncertainty in N, whereas the right hand panel illustrates a specific (and baroque) 
post-inflationary history (solid line) which then flxes the effective expansion history 
and w. 



the (supposedly) post-inflationary evolution. The definition of w can accommodate a 
secondary burst of infiation, but the overall evolution must be such that modes are 
re-entering rather than leaving the horizon. If the universe thermalized instantaneously 
and the number of degrees of freedom in the thermal bath does not change as the 
universe expands, w = 1/3. Conversely, w = implies the universe is (on average) 
matter dominated until p^/^ = pj^f^. 

Note that preh usually refers to the actual energy scale of thermalization, whereas 
in this treatment preh is the energy scale at which thermalization is guaranteed to have 
occurred. As noted above, direct experimental constraints require pj^^ > C(10) MeV. 
Theoretical considerations may be used to justify a much higher value of p^^^, but it is 
important to draw a clear distinction between a theoretical prior and an observational 
bound. As we illustrate in Figure El the value of iD is a function of preh, given a fixed 
cosmological evolution. 

Equation ( IT9|) appears in the literature in various forms. Assuming a matter 
dominated phase followed by thermalization {w = 0) one can derive the result quoted 
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where preh is again the usual thermalization scale. This expression assumes the effective 
number of thermal degrees of freedom is constant as the universe expands, and makes 
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no mention of the contribution of neutrino mass or dark energy to the expansion of the 
universe after nucleosynthesis. 

We distinguish three priors for the evolution of the post-inflationary universe. The 
first, and most restrictive, is to set w = 1/3 and stipulate that the universe is thermalized 
at all scales between the end of inflation and matter-radiation equality. Secondly, we 
can posit a matter dominated phase, followed by thermalization, so w = 0, and preh is a 
free parameter. Finally, if we foreswear all knowledge of the post- inflationary expansion 
history that is not astrophysically verified, w) is a free parameter with — 1/3 < u) < 1. 

It is useful to view AA^ as a "shift" in the pivot /c^, which is mapped to k = k^e^'^ , 
where k is measured in a universe which thermalized immediately after inflation. 
Even with the "middle" assumption of matter domination followed by thermalization, 
AN{k) ~ —9 if {pend/PrehY^^ ~ 10^^, relative to instant reheating. Anticipating the 
results of the following section, simple inflationary models have 10~^ < \as\ ^ 10~^. 
Recalling the definition a^ = dus/dlogk, the corresponding theoretical uncertainty in 
Ug is (conservatively) \6ns\ ~ 5 x 10"'^, without any appeal to exotic dynamics in the 
early universe. For many inflationary models this theoretical ambiguity can easily exceed 
the forecast statistical uncertainty in Ug expected from Planck |31]. Finally, with the 
full range of w, equation (TT9l) gives —25 < AA^ < 8.4 if we insist on thermalization at 
the TeV scale, and an even greater range if we set preh at the MeV scale. 

3. Inflationary Predictions 

We now study the connection between w and the predictions of specific inflationary 
models. Given a candidate model of particle physics which makes a full (and 
computable) set of predictions for the evolution of the universe, w would not be a 
free parameter. In practice, however, inflationary models are specified independently of 
a full theory of particle physics. Consequently, the predictions of inflationary models 
should properly be considered as a combination of the inflationary dynamics and w. 
We begin by collating the predictions for r, Ug and a^ of several broad classes of simple 
two parameter inflationary models. In each case, we eliminate one free parameter from 
the model by matching the amplitude of the perturbation spectrum at the pivot to the 
observed value, Ag in the current concordance cosmology. In practice, the following 
results depend very weakly on the precise value of this parameter. 

3.1. Single term potentials: 0" 

We start with models whose potentials are simple powers of 0, 

V = AM^7'^^, (21) 

n 

so A is a dimensionless constant. Once upon a time, one might have assumed that 

n is an even integer, for which n = 2 is currently the only viable value. However, 

recent developments in string cosmology yield models with n = 2/5, n = 2/3 and 
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n = 1, so we treat n as a continuous quantity [71 |H]. In these models, inflation ends at 
<Pend = nMpi/\/2 and the field value N e-folds before the end of infiation, 07v, is 
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We can write the slow roll parameters as a function of A^ 
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so the spectral parameters are (at lowest order in A^) 
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{AN + nY AN + n {AN + n)^' 

Note that r — )■ as n — > 0, but a^ — ^ —l/N"^ in the same limit. Given that N 
is bounded above, the running remains non-zero in these models, even as the tensor 
component becomes vanishingly small. 



3.2. Natural Inflation 

Natural infiation [6j is governed by the axion-motivated potential. 
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Setting e = 1 we solve for the end-point of infiation. 
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Making no assumptions about f/Mpi one finds 
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When / is small relative to Mpi, infiation ends when cos(0//) is close to unity - in 
other words, at the top of the hill. Conversely, when / is substantially larger than Mpi, 
infiation continues until the field is close to its minimum and cos {(f) end/ f) -^ —1- In 
this case the final 60 e-folds of infiation are generated as approaches the bottom of 
the axion potential: e, rj and ^ are all positive, yielding an observable tensor spectrum 
and a moderate running. Conversely, for smaller /, astrophysically relevant modes are 
generated at smaller values of (j). For extremely small /, e, and ^ ^ \ri\ in this limit, 
and we have a vanishing running, a vanishing tensor amplitude, and ra^ — 1 is fixed by 
rj. The lower bound on n^ sets the minimal value of /: the three year WMAP data set 
gives / > 0.7\/87rMpi or / > 3.5Mpi [IS]. Interpolating between these limits suggests 
that as{ki,) approaches zero as / decreases, and this expectation is confirmed in the 
numerical analysis that follows. 
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Figure 3. Top panel: Ug vs. r, Bottom, Ug vs. as, as vs. r. We plot </>"■ inflation 
(red), natural inflation (green), inflection (purple) and hilltop (blue). Hilltop and 
inflection point inflation meet the (/>" curve at the point where n — 1; natural inflation 
is degenerate with m'^ip'^ inflation in the limit that / is very large. 



3.3. Hilltop and Inflection Point Inflation 



Now consider the potential 



n 



(2^ 



where A is a dimensionless constant. This potential can be specialized to several distinct 
models: e.g. hilltop (n = 2 and n = 4) and inflection point {n = 3). Like natural 
inflation, these models have two distinct limits. In one limit A is relatively small, and 
almost all inflationary growth occurs near = 0, with e <t^ rj. Alternatively, if A^ is 
large (relative to the GUT scale), the physically relevant portion of inflation occurs with 
(f) far from the origin, and e ~ |?7|. In the latter case the potential is effectively linear, 
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Figure 4. Same models as Figure |31 assum.ing a period of matter dominated expansion 
before thermalization, with Trek > 10^ GeV. Colors denote the (logarithmic) reheating 
temperature - orange is low, and blue is high. 



and V{(f)) oc - or 0", with n = 1. If we restrict attention to "small 0" we find 
N ^ ^' ^ 



(n-2)M^f" A 0"-2' 
To lowest order we can compute 



A 



12 



127r2M^f-2n ^2 



{n-2)XN 



(2n-2) 
(n-2) 



(29) 



(30) 



When n = 4, A drops out of the above expression and inflation can occur at any 
energyl5ll Consequently, for n = 4, r is essentially a free parameter. With n = 3, we 

^ To be strictly accurate, the energy scale enters logarithmically via N when we solve the constraints 
self-consistently. 
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Same models as Figures [3] and IH but for the full range of w and 
Pr..h — 1 TeV. Parallel lines for each model denote fixed values of A^. 



Figure 5 

1/4 
reh 



recover the inflection point model ^\. Current data suggests n^ < 1, so we can ignore 
the singularity in the slow roll expressions at the = 0, since red spectra require > 0. 
For these models, as e becomes small, ^ grows large. Consequently, in the low r limit 
we find a substantial, negative running for the hilltop and inflection scenarios. 



3.4- Observables and the Post-Inflationary Thermal History 

The above catalogue made no direct reference to the post-inflationary evolution of the 
universe. Figure [3] shows the parameter values derived assuming instant thermalization 
at the end of inflation. Each model has two free parameters to the potential, one of 
which is removed by matching to the observed amplitude of the perturbations at CMB 
scales. The form of these plots is very weakly sensitive to the precise value we take for 
this parameter, so these plots deflne one-dimensional curves in the parameter space. 
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Figure 6. We now plot the ranges of ris and ag given by our set of models when 
r < 0.05 (left) and r < 0.01 (right). 



If we assume that thermalization is preceded by a matter dominated phase, and that 
reheating has occurred before the universe reaches TeV scale densities, > AA^ > —9, 
relative to the instant reheating case. Figure Hlreplots the models from Figure [3] including 
the effect of the unknown thermalization scale. A nontrivial portion of the parameter 
space is now occupied, and each parameter plane contains degeneracies which are not 
present for instant preheating. Finally, Figure |5] plots the inflationary observables with 
Preh ^^ ^^^ '^^^ scale and arbitrary w. In the absence of a theoretical prior for the 
post- inflationary evolution these simple models are significantly degenerate, and cover 
a large fraction of each parameter plane. 

We can also ask how these plots change when we have more information about some 
of these parameter values. For instance, assuming that Planck and the next generation 
of sub-orbital polarization experiments put a tight upper bound on r: Figure [6] shows 
the permitted regions of the {ns,as) plane for our models that satisfy r < 0.05 and 
r < 0.01. Further, in either of these scenarios Ug is also likely to be tightly constrained, 
to the point that at least one of the classes of models we consider will be ruled out. 

For this collection of models, r varies by orders of magnitude, but typically 
~10 ^ < c^s < — 10~^, so while a^ is small, it is far more consistent. Detecting a 
running at the lower end of this range will be a challenging task as we see in the next 
section, but this range provides a fairly clear target for future observational campaigns. 



4. Observational Constraints 

In this section, we explore how we might probe the running at the 10^^ > la^l > 10~^ 
level. As a base level we estimate the constraints from Planck [H] using the Fisher 
matrix formalism [17] and using instrumental parameters taken from [3]. These 
calculations suggest that Planck should achieve constraints on the tilt and running of 
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Sus = 0.003 and Sas = 0.005 respectively. To improve upon this will require information 
from large scale structure, which may come from galaxy surveys, 21 cm experiments, 
cosmic shear surveys, or the Lyman alpha forest. We will consider only the first two of 
these probes here as having the most potential to provide constraints on models of slow 
roll inflation. 

Cosmic shear surveys will be very powerful in constraining dark energy models, 
but provide only limited improvement to constraints on inflationary parameters |48j . 
Perhaps a factor of two improvement on us and as beyond constraints from Planck 
level data seem possible |l9l [50]. Several attempts have been made to use the Lyman 
alpha forest to constrain inflation using spectra from SDSS [511 |52l |53] . While this could 
be very powerful since the forest constrains scales k ~ l/iMpc~^, greatly extending 
the lever arm, modeling the ionization state and thermodynamic properties of the 
intergalactic medium to convert flux measurements into a density power spectrum is very 
challenging [51]. The Baryon Oscillation Sky Survey (BOSS) is expected to increase the 
power of Lyman alpha forest constraints by greatly increasing the number of quasar lines 
of sight and allowing non-parametric reconstruction of the inflationary power spectrum 
[55] . Robust estimation of inflationary constraints from the Lyman alpha forest is an 
important topic, but it lies beyond the scope of this paper and we leave it to future 
work. 

4.1. Galaxy surveys 

4-1.1. Fisher matrix formalism Cosmological constraints from galaxy clustering and 
baryon acoustic oscillations (BAO) play a crucial role in breaking degeneracies in the 
CMB [561 EZ] • The constraining power of a galaxy survey is determined by three main 
factors: the volume of sky surveyed V^urvey, which determines the number of k modes 
that fit within the survey; the number density of galaxies observed rigai, which determines 
the importance of shot noise; and the maximum wavenumber where the analysis ends 
^max, usually set by the scale on which structure becomes non-hnear. Since the non- 
linear scale is a strong function of redshift, we must also consider the redshift range 
(zjnin, -2max) of the survcy. It has been shown [58] that high redshift galaxy surveys can 
reasonably constrain the running at the level as ~ 10^^. We explore both this scenario 
and more ambitious hypothetical examples to probe the limits of the possible. 

Following the basic formalism set forward in [56], we calculate the constraints for 
Planck in combination with galaxy surveys. The Fisher matrix for a galaxy survey can 
be written as [56] 



, d^k (91ogP(k)91ogP(k) 



"^ ' 2(27r)3 dpi dpj 



%alP(fc, ^) 1 ^ y^^^ ^ .g^. 

'0 ^\^n J upi ujjj \_ngs.iP{k, fl) + 1\ 

where the derivatives are evaluated using the cosmological parameters of the fiducial 
model. The 1 — a errors on the parameter pi are then given by Api = \l F~^^ . 

In calculating the power spectrum of galaxies, we allow for biasing and redshift 
space distortions, so that P{k,fi) = (1 + P^'^^Pgik) + Pq, where /i = k\\/k, Pq is 
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residual shot noise, and the distortion parameter /3 = Qm{z)^'^/bi- This gives the 
intrinsic power spectrum, and the k space power spectrum inferred from angular and 
redshift measurements is further modified by the Alcock-Paczynski effect [59]. We 
assume that on large scales the galaxies are linearly biased with respect to the dark 
matter distribution so that Pg{k) = h\Ps{k), with linear bias hi. Nonlinear biasing can 
introduce considerable systematic uncertainty into the galaxy power spectrum on small 
scales, and we discuss this in more detail below. 

In addition to the non-linearity of galaxy bias, we should worry about the non-linear 
evolution of the underlying density field which can erase cosmological information on 
small scales. We initially model the cut-off scale fcmax using the prescription of [60] , who 
assume that non-linear effects are important for k > kni = Tc/{2Rni), where Rni is the 
scale on which averaged density fluctuations a{Rni) = 0.5. Beyond /cmax we are throwing 
away cosmological information, which a better understanding of non-linear effects (e.g. 
via higher order perturbation theory) might allow us to exploit. Consequently, we 
generalize this definition to kni{'y), where a{Rni) = 7, and use 7 = 0.5 unless otherwise 
specified. 

In calculating the Fisher matrix, we include the cosmological parameters [Vtmh?, 
flhh'^, Q\, As, T, Us, as, r). To this we add parameters describing the galaxy bias, hi, 
redshift-space distortions, /3, and shot noise, Pq, in each redshift bin. The fiducial value 
for the bias is calculated using the Sheth-Torman peak-background split [61] appropriate 
for dark matter halos and assuming that each halo hosts only a single galaxy (see [62] 
for a review of the technical details). We set the minimum halo mass Mmin using the 
mean number density of galaxies in the survey rzgai by assuming that the number density 
of halos rahaio(^ > ^min) = 'T-gai- This is appropriate, for example, if the survey collects 
only the brightest object within the survey volume. 

Our numerical calculations and forecasts use the fiducial cosmology Vtm = 0.3, 
^A = 0.7, VLb = 0.046, H = lOO/ikms^^ Mpc^^ (with h = 0.7), ns = 0.95, r = 0.1, and 
o"8 = 0.9, consistent with the WMAP seven year data [63] • We note that our fiducial 
set of cosmological parameters is not exhaustive and could be expanded, for example to 
include the dark-energy equation of state w and the sum of neutrino masses M^, as free 
parameters. Including more parameters will degrade the constraints somewhat beyond 
those considered here. 

4- 1-2. Constraints from galaxy surveys We first delineate the parameter space 
accessible to galaxy surveys to get a sense of how well present and future galaxy surveys 
might perform. We initially focus on the Q — Ugsi plane, where Q is the observing area, 
noting that there is a tension between going deep, to increase ngai, or wide, to increase fl. 
As a point of reference, the HUDF [Hubble Ultra Deep Field], which covers only 11 sq. 
arcmin., has a number density of galaxies rigai = 0.02h^Mpc~^ at z = 5 [61], obtained 
from integrating the fitted luminosity function of drop outs down to the observational 
magnitude limit of M = —16. This gives a sense of the upper limit for the number 
density of galaxies accessible to a very futuristic survey. To a similar brightness limit. 
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Figure 7. Contour plot of Arts (left panel) and Aag (right panel) at the 1-cr level in 
the $7 — Tigai plane at z = 3 for a galaxy survey in combination with Planck. 



about twice this number of Lyman break galaxies are observed at 2 = 3 [65] . 

We show the constraints on Us and a^ in Figure [3 for different values of (fi, rigai) 
for a galaxy survey centered at 2; = 3 of width A2; = 1, covering an angular region 
VL. Increasing Vt always improves the constraints, since more volume gives greater 
statistical information. Increasing ngai also improves the constraints, but only up to 
the point that shot-noise becomes negligible. Galaxy surveys have a larger impact on 
as than n^, since Us is already tightly constrained from CMB data alone. To give 
a sense of what future surveys might achieve. Figure [7] includes points denoting the 
expected performance of several upcoming or proposed spectroscopic galaxy surveys 
that will reach 2 = 3, including WFMOS with (fi,ngai) = (300, lO'^) (in units of 
(arcmin^, /i^Mpc"^)), and HETDEX with (fi,ngai) = (420,0.5 x lO'^). Moreover we 
present results for two possible configurations of the Cosmic Inflation Probe (CIP)l1j, a 
proposed satellite mission designed to observe Ha galaxies over the range z = 1.9 — 6.6, 
CIP3 with (n,ngai) = (10^,7 x 10~^) and a somewhat beefed up version, CIP4 with 
(n,ngai) = (10^,1.5 X 10^'^). CIP is designed to reach higher in redshift to ;2 < 6.5, 
which would further improve its inflationary constraints. 

Extra volume can be added by expanding the redshift range covered by the survey. 
In Figure [H we show the parameter constraints for an all sky galaxy survey in which 
we add redshift bins of width Az = 1 up to a maximum redshift -Zmax- Observations 
at higher redshift benefit from probing both more volume (for fixed Q), and because 
shorter comoving scales are still in the linear regime, leading to tighter bounds on n^ 
and ag. However, the extra linear scales are only useful if shot noise is not significant, as 
can be seen in the saturation of the curves in Figure [8] for fixed ngai- In this calculation, 
the constraint is usually dominated by the highest redshift bin. 

It is apparent that galaxy surveys focussed on redshifts z < 3, such as SKA or 



http://www.cfa.harvard.edu/cip/ 
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Figure 8. Constraints on inflationary parameters for Planck in combination with a 
full sky survey up to redshift z,„ax for "-gal = 10^* — 10^^ h^Mpc^^ . Black dotted 
curve is the cosmic variance limit. 



JDEM, will not detect running at levels below la^l ~ 10^^. Pushing into the range 
Aag = few X 10~^, characteristic of the slow roll models considered here requires reaching 
higher redshifts. Although these might be accessible to a space based survey, attaining 
such redshifts may require different techniques, such as 21 cm experiments. 

Alternatively, rather than pushing to higher redshifts we might hope to make better 
use of data at scales in the non-linear regime. Surveys such as SDSS already measure the 
power spectrum on non-linear scales, but are unable to make full use of this information. 
In principle, quasi-linear modes may be used to obtain cosmological information if they 
can be accurately modelled by, for example, higher order perturbation theory [661 EZ] ■ 
In Figure [HI we show the improvement in a survey that reaches to 2; = 3 if modes out to 
^max(7) can be included. The Ug constraint is relatively insensitive to k^^x, while gains 
in as by a factor of a few can be achieved by accessing smaller scales. This agrees with 
the findings of ^8] . As expected, the flattening of the contours with increasing 7 shows 
that to take advantage of the extra scales larger rigai is required. 

While an understanding of nonlinearities may improve our ability to use small 
scale modes, non-linear or non-local biasing may prevent this. Non-linear biasing 
may be modeled via the halo occupation distribution (HOD) [6H] or via higher 
order perturbation theory [69j. However, either approach introduces additional free 
parameters that must be marginalized over, degrading parameter constraints. Inclusion 
of higher order bias terms has been shown to have relatively little effect on baryon 
acoustic oscillations [70j, where the oscillatory features are robust to changes in the broad 
band power, but may prove more important when constraining inflationary parameters. 
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Figure 9. Contour plot of Aug (left panel) and Aa^ (right panel) in the 7 — rigai 
plane at z = 3 for il = 10, 000 sq. deg.. 



At the level of the Fisher matrix, we have found that accounting for third order bias 
parameters, of the form detailed in [SH], degrades the constraints on Ug and a^ by a 
factor of ~ 2. Systematic effects from the non-linear biasing are beyond the scope of 
our analysis. Additionally, non-local biasing can occur if, for example, the large scale 
features of reionization modulate the galaxy power spectrum via suppression of galaxies 
in ionized regions [71]. This would also modify the broad band power of the galaxy 
power spectrum introducing systematic uncertainties. 



4-2. 21 cm Experiments 

In the preceding section, we have shown the need to utilise high-redshift galaxy surveys 
in order to push constraints on Ug to the 10"'' level. These redshifts may also be 
probed with the redshifted 21 cm signal which, in the absence of ionization fluctuations, 
tracks the density field. Several low-frequency radio interferometers including GMRT 
[72], LOFAR0, MWA0, and PAPER [73] are currently under construction and hope 
to make an initial detection of this signal. While these first generation experiments 
are unlikely to add much to our knowledge of cosmological parameters, future 21 cm 
experiments such as SKA lffl have the potential to significantly improve upon our current 
knowledge of inflationary parameters. The analysis of [71] showed that a futuristic 21 
cm experiment, the Fast Fourier Transform Telescope (FFTT) [75], could potentially 
achieve Aa^ = 10~^, under optimistic assumptions about reionization. 

Such experiments can be focussed either at moderate redshifts z = 7 — 25, where 
the effects of reionization and the first stars will serve as a major contaminant [76j, or 
at redshifts z = 25 — 50, before star formation has begun. Observing high redshifts 

* http://www.lofar.org/ 

tt http://www.MWAtelescope.org/ 

ft http://www.skatelescope.org/ 
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requires removing foregrounds, especially galactic synchrotron emission which typically 
scales as z/~^'® and dominates the signal. We consider both regimes separately. We 
also note that 21 cm intensity mapping fTll 178] potentially offers another way of going 
after the density power spectrum. Since the challenges and redshifts associated with 
intensity mapping are largely the same as those associated with galaxy surveys, we will 
not consider them in detail. 

4-2.1. Cosmology from the Epoch of Reionization Observations of the 21 cm signal 
from the epoch of reionization (EoR) constrain a combination of the fluctuations in 
the density field and the ionization field. Neglecting fluctuations in the 21 cm spin 
temperature (which is a reasonable, although not guaranteed to be safe [76j), the 21 cm 
brightness temperature power spectrum takes the form 

Pt, = Pss + 2Pfe + P.. - 2fi^ {Pss + P.s)+ li^Pss + Pf(,k,^A ■ (32) 

where k and the line of sight, Pss, Pxx, and Pxs are the power spectrum of the density 
field, the ionization field, and the density-ionization cross-correlation respectively, while 
fi is again the angle between the Fourier mode and the line of sight. The final term 
Pf(k,fjL) contains contributions higher than quadratic in the perturbations. Although 
these would normally be neglected as small, during reionization fluctuations in the 
ionized fraction Xi can be of order unity, so terms higher than quadratic contribute to 
the power spectrum. These terms spoil the simple angular dependence expected from 
linear theory [79] and potentially degrades our ability to separate astrophysics from 
cosmology beyond that considered in [71]. Numerical simulations [801 E] show that 
higher order terms contribute significantly to the power spectrum on all scales once 
reionization is underway, and must therefore be included if cosmological parameters are 
to be correctly estimated. 

We model Pxx and Psx using fitting functions and fiducial parameters from |74j . 
These take the form 

Psx = bl^exp[-axs{kRxs) - {kRxsfjPss 

Pxx = &L[1 + oixx{kRxx) + [kRxxf]-"-^ Pss, (33) 

where hsx, O-xs, RxS, bxx, o^xx, Rxx and '-)xx are free parameters to be varied in our Fisher 
analysis. 

In our Fisher analysis, we will consider two scenarios for the ionized fiuctuations: 
an optimistic scenario where ionization contributions are neglected (OPT), and a case 
where only quadratic contributions to the power spectrum are included (MID). We have 
separately considered the more detailed case in which higher order terms are considered. 
These require considerably more work and appear to degrade the inflationary constraints 
by only an additional ~ 50% over the MID case. For this reason, we leave a more detailed 
discussion of the higher order terms to separate work. 

4.2.2. Fisher matrix formalism We perform a Fisher analysis of 21 cm experiments 
following [82l [79]. The Fisher matrix for a 21 cm experiment is given by a sum over 
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at z = 8 in the OPT scenario. 
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where the variance of a 21 cm power spectrum estimate for a single k-mode with line of 
sight component k\\ = fik is given by 
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This depends upon the system temperature Tgys, the survey bandwidth B, the total 
observing time tint, the conformal distance D{z) to the center of the survey at redshift z, 
the depth of the survey AD, the observed wavelength A, and the effective collecting area 
of each antennae tile A^,. The effect of the configuration of the antennae is encoded in 
the number density of baselines n{k±) that observe a mode with transverse wavenumber 
k±. Observing a number of fields A^fiew further reduces the variance. We will consider 
only the case of a filled array and note that for baselines well below the cutoff due to the 
finite size of the array n{kj_) fa A^antA^/Ag = (AtotMe)(A^/Ae). The response of a single 
dipole sets the minimum effective area A^ = A^/4. These arrays are therefore usefully 
described in terms of Ae, Atot, Bti^t, and the redshift range covered. 

The sensitivity of 21 cm experiments in the Ae — Atot plane for i? = 8 MHz, 
Afield = 2, and 4000 hours integration time is shown in Figure [101 ^tot affects both 
the sensitivity and the angular resolution, which controls the largest kj_ accessible. 
Ae also impacts the sensitivity and controls the volume probed by the survey via the 
instantaneous field of view. Once the experiment becomes sample variance limited the 
only way to improve sensitivity is via increasing the volume probed, i.e. by making Ae 
smaller. 
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Figure 11. Contour plot of Aug (left panel) Aog (right panel) in the A^ — Atot plane 
at z = 8 in the MID scenario. 



It was shown in [71] that the two most important hmiting factors for 21 cm 
observations are, respectively, modehng of the contribution to the power spectrum 
from ionization fluctuations and foregrounds. In Figure [TOl we show contours in the 
case where ionization fluctuations are ignored (OPT), while in Figure [TT] we allow for 
quadratic terms (MID). In the MID case constraints are degraded by a factor of ~ 2 
over the OPT case. 

In each of these figures, we mark the approximate position of several current and 
future 21 cm experiments. We indicate MWA, built out of 500 antennae tiles with 
Aes = 14m^ at 2; = 8, and two possible successor instruments built by increasing the 
number of tiles by a factor of 10 (MWA5000) and 100 (MWASOk). We also mark the 
specifications of the proposed FFTT, composed of 10^ dipoles with A^s = A^/4. We 
have assumed a compact array design for all of these, which is not actually how MWA 
is being built, although this give a small correction. MWASOk and FFTT have similar 
collecting area, the boost from cross-correlating all of the dipoles makes FFTT much 
more sensitive to the power spectrum, both from a greater raw sensitivity and since it 
surveys a larger volume of sky. 

As expected from our earlier discussion, the increased volume accessible from a 
full sky survey at 2; = 8 allows constraints on the inflationary parameters at the level 
needed to detect as at levels which are characteristic of simple slow roll models, and 
FFTT has the instrumental sensitivity required to actually make measurements at this 
level. Allowing for ionization contributions degrades the sensitivity to the running, but 
the high sensitivity of the instrument means that even after this degradation, interesting 
constraints on the running are still obtainable. However, the need to model ionization 
contributions raises the possibility of systematic biases, especially at this high level of 
precision and it is still unknown whether future 21 cm experiments will be able to remove 
foregrounds and control systematics at the level needed to achieve this sensitivity. 
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Figure 12. Contour plot of Aug (left panel) Aag (right panel) in the Ac — Atot plane 
at z = 30 in combination with EPIC. 



4-2.3. 21 cm cosmology from the dark ages We have seen that trying to obtain 
cosmology from the 21 cm signal during the epoch of reionization is complicated by 
the presence of fluctuations in the neutral fraction. In principle, one can avoid this by 
going to redshifts 2; > 30, before star formation has begun. Here brightness fluctuations 
are expected to trace the density as a result of collisional coupling [831 EH] • 

There are two main challenges to accessing this redshift range. The first is the 
Earth's ionosphere, which has a plasma frequency of ~ 10 MHz leading to smearing 
of the signal and a resulting loss of angular resolution that prevents observations of 
the z > 30 regime from the ground |8l]. This requires that observations be conducted 
from space, including proposals based on the lunar surface. The second challenge is 
that foregrounds scale with frequency as v~'^-^ and become very large at the frequencies 
of interest, making observations significantly more difficult than during the epoch of 
reionization [76]. Proposals for lunar arrays exist, such as LARC and DALI [85] , 
although these are geared towards arrays with size ~km^. Since low frequency dipole 
antennae can be little more than wires printed on a sheet of plastic a lunar array is 
less crazy that it might initially appear, and the proposed Ares V heavy launch vehicle 
could deliver ~ 0.5 km^ of collecting area to the lunar surface. 

In Figure [121 we explore the requirements of measuring the tilt and running with 
a lunar array. Given the time scale required for building such an array, we consider the 
inflationary constraints from a 21 cm instrument in combination with EPIC [86j, which 
by itself achieves An^ = 0.0018 and Aa^ = 0.0026, in our calculations. Arrays with 
collecting area of order 10 km^ are required to measure the 21 cm signal and achieving 
a high signal of noise requires even larger collecting areas. An FFTT like array with 
collecting area 10^ km^ could measure the running at the level of 10~^ from the lunar 
surface, which is clearly an extremely futuristic proposal. Smaller arrays might still 
be useful for constraining inflation via observations of the large scale power spectrum, 
which can be used to constrain compensated isocurvature modes [87] . 

We list the constraints on the inflationary sector that can be obtained with our 
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Table 1. Inflationary parameter constraints for CIP and FFTT. The first block gives 
forecasts for CIP and FFTT on their own, the second block gives forecasts for Planck, 
Planck+CIP and Planck + FFTT. The Fisher forecasts are performed with the fiducial 
parameters Ug = 0.95, as = —0.0005, and r = .0045. 



Experiment 


An, 


Aas 


Ar 


CIP3 


0.0089 


0.0028 


- 


CIP4 


0.0028 


0.0009 


- 


^^__ OPT 
FFTT 

MID 


0.0011 


0.00023 


- 


0.00082 


0.00032 


- 


Planck 


0.0032 


0.005 


0.058 


+CIP3 


0.0019 


0.0011 


0.05 


+CIP4 


0.0011 


0.0006 


0.048 



^^^^ OPT 0.00034 0.000095 0.048 
+FFTT 

MID 0.00067 0.00028 0.048 



fiducial versions of CIP and FFTT in Table [H CIP3 (CIP4) assumes a galaxy survey 
over 1000 (10000) sq. deg. with redshift bins at z = 3, 4, and 5 achieving galaxy 
densities of rigai = 8 x 10~^, 4 x 10~^, and 1 x 10~^ h^Mpc^ respectively. For FFTT, we 
assume 10^ dipoles with A^^ = A^/4, for a collecting area of Atot ~1 km^ aX z = 8. For 
this illustration, we take the fiducial inflationary parameters corresponding to natural 
inflation with / = a/Stt and A^ = 51. 

In principle, these experiments are capable of precise measurements of the tilt and 
borderline detections of the running expected for simple slow roll models. However, 
once systematic effects from biasing or reionization are taken into account it is apparent 
that a precision measurement of the running is an extremely challenging project. The 
best hope for improving this picture is probably to correctly model modes within the 
non-linear regime, allowing them to included in the analysis. Nonetheless, it is clear 
that there is a sufficient number of comoving modes within the visible Universe that a 
detection of the inflationary running is feasible. 

5. Constraining the post inflationary universe 

The previous section assessed the ability of galaxy surveys and 21 cm experiments 
to detect the running predicted by canonical single field inflationary models. If 
the primordial perturbations were generated by one of these models, the analysis of 
the previous section shows that a direct measurement of a, is possible with future 
observations, given exquisite control of foregrounds and systematics. Figure [13] shows 
the error forecasts for measurements of Ug, r and a,, with a "modest" (Planck- level) 
constraint on the tensor signal, and tight constraints on Ug and a,. For this scenario, 0"' 
models are excluded at high significance. However, while FFTT is forecast to provide a 
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Figure 13. 68-% confidence ellipses for Planck+CIP (dashed ellipse) and 
Planck+FFTT (solid ellipse) drawn on the inflationary parameter space. The Fisher 
analysis is performed assuming natural inflation, and the (rig, as) plot only includes 
models with r < 0.12. We see that in this scenario, all <j>" models are eliminated at 
high significance, but that we cannot use constraints on q^ to break the degeneracy 
between the remaining models. 



marginal detection of as, it cannot break the degeneracy between the remaining hilltop 
models and natural inflation. Conversely, a post-Planck polarization mission would 
detect the tensor signal from natural inflation, and exclude the hilltop models. 

However, as discussed in Sections 2 and 3, the running makes its presence felt 
at much lower levels of precision by inducing a correlation between Ug and the post- 
inflationary expansion history. Consequently, given an explicit inflationary prior, 
accurate constraints on the power spectrum via large scale structure surveys provide 
tight constraints on N and the four variables needed to describe the primordial 
perturbations are typically functions of a smaller number of free parameters in the 
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Table 2. Constraints on reheating temperature and inflationary model parameters 
from Planck, CIP3 and FFTT (MID) for both natural and </>" inflation. Rows with "-" 
give error forecasts for the complementary parameter, with the "-" parameter assumed 
to be specified in the prior. In each case the final row gives a constraint assuming that 
r is measured to ±0.01. 





Natural 


0" 




A^ 


/ 


A^ 


n 


fiducial values 


51 


v^ 


51 


2 


Planck 


5.1 


- 


3.6 


- 




- 


0.33 


- 


0.25 




14.5 


0.93 


19.7 


1.4 


+ ar = 0.01 


3.5 


0.26 


8.6 


0.41 


CIP+Planck 


1.69 


- 


1.2 


- 




- 


0.11 


- 


0.09 




13.7 


0.87 


14.5 


1.14 


+ ar = 0.01 


2.8 


0.18 


3.96 


0.27 


FFTT+Planck 


0.41 


- 


0.29 


- 




- 


0.027 


- 


0.024 




7.0 


0.45 


11.0 


0.91 


+ ar = 0.01 


2.5 


0.17 


2.95 


0.24 



inflationary potential. The reduction in the dimensionality of the parameter space 
imposed by an inflationary prior thus allows us to constrain A^, even when we cannot 
measure a^ with confidence. Written in terms of the inflationary parameters, the Fisher 
matrix becomes 

dqi ^ dqj 






dpn dp, 



(36) 



where the qi are the spectral parameters Ug, a.^, and r while the p„ are the free parameters 
in the speciflc inflationary model, including the number of e-foldings A^. Table |2]lists the 
constraints on inflationary model parameters and A^ that can be expected from Planck, 
CIP and FFTT (in the MID scenario). In practice, these forecasts are approximate 
since the inflationary prior puts sharp cuts on the parameter space, while the Fisher 
matrix assumes a Gaussian likelihood. 

Given a tight prior on the inflationary potential, we see that A^ is likely to be well 
constrained by Planck, so that Planck could begin to distinguish between quadratic 
inflation followed quickly by an unbroken period of radiation domination, and a scenario 
that involved thermal inflation. This knowledge is of signiflcant value to particle 
theorists and inflationary model builders, since it implies a non-trivial correlation 
between the assumed form of the inflationary potential and the properties of the 
post-inflationary universe. Beyond Planck, FFTT and CIP would provide exquisite 
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Figure 14. Top: 68-% confidence ellipses for Planck+CIP (dashed ellipse) and 
Planck -l-FFTT (solid ellipse) for 0" inflation mapped into the N — p plane. Bottom: 
As above, but assuming r is measured to within 0.01. 



constraints on the post-inflationary universe, for a given inflationary scenario. 

Clearly, if the inflaton potential has more free parameters, constraints on A^ weaken 
substantially. For the examples here, estimating both f or n (for natural and (f)"" 
inflation, respectively) and A^ from data leads to substantial play in the allowed values. 
However, a strong constraint on r and a ClP-class measurement of large scale structure 
is sufficient to put a tight constraint on both the shape of the potential and A^. Figure 
[H] shows the constraint forecasts for n and A^, for (j)"' inflation. There is considerable 
degeneracy between these two parameters that neither CIP or FFTT are able to break 
on their own. However, a precision measurement of r breaks this degeneracy, so that 
the parameters are accurately determined by the data. 
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6. Discussion 

We have surveyed the impact of a running spectral index on the predictions of 
inflationary models. Beyond the obvious role of «« as a free parameter in its own right, a 
running index and unknown expansion history of the post-inflationary universe induce a 
substantial theoretical uncertainty in the predictions of inflationary models. Physically, 
the field value at which perturbations are generated is sensitive to the post inflationary 
expansion history, which determines that rate at which modes reenter the horizon after 
inflation. Thus if the universe does not immediately thermalize at the end of inflation, 
the pivot scale is effectively shifted (relative to that defined for instantaneous reheating) 
to ki = he^^. 

If Us was actually constant (as is the case for power-law inflation) changing /c^ would 
not change n^. However, even simple inflationary models have a.^ large enough to ensure 
that the post-inflationary expansion history has a nontrivial impact on the measured 
value of the spectral index. Making the conservative assumption that inflation is followed 
by an effective matter dominated phase with thermalization occurring at the TeV scale 
or above, the resulting uncertainty in n^ is comparable to the statistical uncertainty 
expected from Planck. Conversely, with only minimal theoretical assumptions about 
the properties of the post-inflationary universe the corresponding uncertainty in Ug is 
similar to the statistical uncertainty found in WMAP constraints. 

This ambiguity arises because inflationary models are typically specified solely in 
terms of their potential, rather than being embedded in a larger theory of particle 
physics. Consequently, positing an inflationary model does not specify the properties of 
particle physics from TeV to GUT scales, or the couplings between the inflaton and the 
"rest of particle physics" . With this information, the post-inflationary expansion history 
of the universe could be calculated, and the uncertainty in the inflationary spectrum 
removed. 

Our analysis shows that futuristic large scale structure or 21cm measurements will 
detect as at the level predicted by simple inflationary models, confirming previous work 
in this area. However, a measurement that is accurate enough to distinguish between 
the values of ag predicted by these models will be extremely challenging. The ultimate 
limits to these observations will likely be set by systematic effects, so that we have dwelt 
only upon best possible statistical limits, leaving the issue of systematics aside. 

Conversely, the uncertainty in Us induced by ag will be important to all analyses 
of the inflationary parameter space following Planck, even if ag is not directly 
detected. This apparent paradox is resolved when we recall that measuring ag directly 
requires estimating at least four independent spectral parameters (in addition to other 
concordance variables), namely Ag, r, and Ug, in addition to ag. Conversely, explicit 
inflationary models typically have fewer free parameters, effectively correlating one or 
more of the otherwise independent spectral parameters. Consequently, given a specific 
inflationary potential, we can constrain the integrated expansion history of the post- 
inflationary universe, as described in Section 5. 
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Intriguingly, many supersymmetric scenarios predict that the primordial universe 
undergoes a long matter dominated phase, due to the presence of "light" moduli whose 
energy density scales like non-relativistic matter. In some scenarios these moduli are 
erased by thermal inflation [38] ; in others they decay of their own accord in such a way as 
to avoid disrupting nucleosynthesis [HE]- In both cases, the number of e- folds required to 
match the pivot scale to the moment it left the horizon during inflation is substantially 
different from that of a universe which is thermalized throughout the post-inflationary 
epoch. Consequently, in the near future astrophysical observations will determine which 
inflationary scenarios are compatible with these common supersymmetric scenarios, and 
which are not. 

The Large Hadron Collider promises to extend our understanding of particle physics 
to the TeV scale, while direct detection experiments are putting significant constraints 
on many dark matter models. If these experiments reveal the mechanisms responsible 
for setting the present-day baryon and dark matter fractions, we will gain significant 
new windows into the very early universe. Given the huge range of energies that can 
lie between the inflationary and TeV scales, along with the lack of direct probes of 
particle physics at these energies, the post-inflationary era amounts to a fundamental 
"dark age" , where both the cosmological evolution and fundamental laws of physics are 
unknown. We have seen here that the relationship between n^ and the post-inflationary 
expansion history induced by a^ may help to illuminate this currently mysterious epoch. 
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